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ABSTRACT 

A  Hi|^-Temperature  Test  Facility  Is  being  installed  in  the  Air  Force  Aero  Propulsion 
Laboratory  for  analjraing  various  sophisticated  cooling  designs  over  an  extended  gas 
temperature  range.  The  objective  is  to  evaluate  these  schemes  with  respect  to  gas  temperature 
and  pressure  levels.  The  experimental  data  will  be  correlated  with  theoretical  predictions  to 
provide  a  basis  for  future  hi^temperature  turbine  blade  designs. 

This  report  describes  a  theoretical  prediction  of  the  convection  heat-transfer  coefficient 
distribution  on  a  test  blade  and  the  computer  program  for  making  the  necessary  calculations. 
The  program  is  written  in  Fortran  IV  ready  for  an  IBM  7094.  The  input  and  output  are 
described.  The  program  starts  with  the  blade  surface  pressure  distribution  obtained 
experimentally  and  yields  the  heat-transfer  coefBcients  on  the  Made  contour. 
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SYMBOLS 


A  channel  cross-sectional  area 

A*  channel  cross-sectional  area  at  which  Mach  number  is  unity 

Cf  turbulent  coefficient  of  skin  friction 

Cp  specific  beat  at  constant  pressure 

h  convection  heat-transfer  coefficient 

H  velocity  shape  factor  (turbulent  boundary  layer) 

n  distance  along  velocity  potential  line 

local  Nusselt  number 

P  pressure 

p_  Prandtl  number 

r 

r  radius 

R  gas  constant 

T  temperature 

V  or  U  velocity 

Vgjy  mid  channel  reference  velocity  at  a  distance  my  from  centerline  of  cascade 

X  surface  distance  from  stagnation  point 

y  distance  along  radial  potential  line 

p  density 

X  ratio  of  specific  heats 

V  kinematic  viscosity 

/X  viscosity 

8  boundary  layer  thickness 

o  displacement  thickness 

o  boundary  layer  energy  thickness 

X  velocity  shape  factor  (laminar  boundary  layer) 

B  momentum  thickness 

r  shearing  stress 
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SUBSCRIPTS 

m  refers  to  tnld-ohannel  oondltlcNis 

CR  refers  to  orltlosl  oomlltion 

o  stage  entrance  conditions  (cold  flow) 

1  stage  entrance  coodltlcNis  (hot  flow) 

n  distance  along  velocity  potential  line 

no  distance  along  potential  line  (cold  flow) 

nl  distance  along  potential  line  (hot  flow) 

S  static  coodltlona 

ad  adiabatic  conditions 
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SECTION  I 
INTRODUCTION 


Currently,  the  Air  Force  Aero  Propulaion  Laboratory  ie  Inatalling  an  in-bouae  High- 
Temperature  Teet  Facility  which  will  be  ueed  initially  to  analyze  varloua  cooling  techniques, 
aa  applied  to  turbine  engine  bladee,  over  an  extended  gas  temperature  range.  The  objective 
Is  to  evaluate  theae  cooling  eohemea  with  reepeot  to  gas  temperature  and  pressure  levels. 
Experimental  data  will  be  correlated  with  the  various  theoretical  predictions  of  blade  cooling 
effectiveness,  thus  providing  a  basis  for  ftiture  hi^-temperature  turbine  blade  designs 

Ibe  experimental  data  used  is  that  obtained  from  the  recording  of  surface  metal 
temperatures  by  Infrared  thermometry.  In  the  theoretloal  prediction  of  turbine  blade 
temperatures,  cooled  or  uncooled,  the  major  steps  are  the  calculation  of  the  convection 
heat-transfer  coefflcients  on  the  blade  eurfhce  and  the  calculation  of  the  driving  or  adiabatio 
wall  temperatures.  Ilieee  calculationo  depend  primarily  on  the  accuracy  with  which  the 
behavior  of  the  boundary  layer  Is  predicted.  It  is  generally  assumed  that  the  Introduction  of 
a  cooling  fluid  In  the  boundary  layer  has  a  negligible  effect  on  the  convection  heat-transfer 
coefflcients,  but  It  affects  considerably  the  adiabatic  or  driving  wall  temperature. 

This  report  presents  a  theoretical  approach  to  the  problem  of  determining  the  adiabatio 
wall  temperatures  and  the  convection  heat-transfer  coefficients  . 

This  program  le  written  in  Fortran  IV  ready  for  use  on  an  IBM  7094  for  an  uncooled 
blade  and  the  computer  program  to  which  it  led.  The  results  obtained  from  the  program 
developed  here  are  compared  with  those  from  the  computer  program  developed  by  the 
Allison  Division  of  General  Motors  on  Contract  AF  33(6 15)- 298 5. 
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SECTION  n 

THEORY  AND  DISCUSSION 


llie  veloolty  of  a  perfect*  frlotlonleee*  oompreeslble  gne  flowing  through  a  channel*  at 
any  point  across  a  potential  line*  ie  given  by  the  equation 

Cm)]  I" 

where*  for  a  nourotatlng  channel* 

Vfny  *  V  •«P  [  /^  "  ^  *2) 

^  'yo 

is  the  mid  channel  velocity  at  a  distance  y  from  a  reference  point  y^  along  a  radial  potential 
line  and  '^yo  Is  the  velocity  at  the  reference  point  (Figure  1).  It  Is  seen  from  Equations  1  and 
2  that  V  Is  strictly  a  function  of  channel  geometry  and  Vyo‘ 

In  the  computer  program  currently  In  use  at  AFAPL  (Reference  1)*  Vyo  is  arbitrarily 
chosen;  then*  the  fluid  velocity  and  density  are  evaluated  at  specified  locations  along  a 
veloolty  potential  line*  the  latter  being  obtained  from  the  equation 


^0 


(3) 


Integration  of  the  product  over  a  given  channel  cross  section  yields  the  total  mass 

flow  which  Is  compared  with  the  required  flow  (usually  the  maximum  or  choked  flow).  Vyg  is 

Improved  until  dA  (required  flow)  <  tolerance.  Use  of  this  computer  program 


requires  an  accurate  description  of  the  channel  configuration  at  each  of  the  cross  sections 
where  the  channel  wall  velocity  must  be  obtained.  Eleven  cross  sections  are  presently  used 
to  provide  enough  data  for  Interpolation  of  the  velocity  profile  along  the  channel  walls 
between  the  stage  entrance  and  the  throat*  beyond  which  extrapolation  is  depended  upon  to 
obtain  the  velocity  profile.  As  inputs  to  the  program*  curvatures*  length  of  velocity  potential 
lines,  and  the  angle  between  the  gas  stream  direction  upstream  of  stage  entrance  and  the  mid 
channel  stream-line  must  be  measured  or  calculated.  A  table  of  V/V^^.  versus  axial  chord 


2 


AFAPL-TR-68-143 


distance  must  be  prepared  from  a  curve  drawn  through  the  extrapolated  and  interpolated 

points  as  an  input  to  a  second  program  (described  later)  which  in  turn  translates  V/V  versus 

CR 

axial  chord  distance  into  versus  surface  distance  from  the  blade  stagnation  point  and 

ultimately  yields  the  convection  heat-transfer  coefficients  along  the  blade  profile. 

It  follows  that,  even  when  only  changes  In  the  gas  flow  parameters  (P^,  T^,  Rq, 
are  involved,  preparation  of  a  new  table  from  the  output  of  the  first  computer  program 
requires  several  hours.  When  a  different  blade  lumflle  must  be  tested,  a  tedious,  time- 
consuming  determination  of  the  channel  parameters  listed  above  must  also  be  conducted. 


A  method  is  herein  described  which  uses  solely  as  Inputs  the  static  pressure  dlstiibutlon 
on  the  blade  surface,  obtained  directly  from  the  cold  flow  test  described  in  Reference  2,  and 
the  gas  flow  parameters  at  t 
method  proceeds  as  follows: 


stage  entrance 

for  both  the  cold  flow  and  hot  flow  tests.  This 

per 

unit 

oreo 

(cold  flow): 

P  f 

,  ()&•»•  1) 

A* 

Rq  To  . 

1 

per 

unit 

oreo 

(hot  flow)! 

1  " 

P  1 

^g  >;  1 

T  r  2  1 

A* 

“  1 

.  R,  T,  J 

1  Xe-iJ 

Fractional  change  In  flow  rate: 


FR  X 


^1  -  Wq 


(5) 


Densities  (inlet  conditions): 


Critical  velocities: 


i  P,  *  \ 


'fcRO  * 


XRI 


2R,  g 


(6) 
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Local  channel  cold  flow  denaittoa: 


p  - 

P 

f'- 

>9-' 

'n 

0 

L 

\ 

'  '  ''CRO  '  J 

Desired  local  channel*  hot-gas 

mass  flow: 

p 

V 

4- 

(P 

V  WfrN 

S  GO 

*00 

no 

noM  } 

Local  channel  hot  flow  density: 

f. 

1-  (. 

— U- 

^  Yno 

■)'l^ 

r 

'  V 

+  1/  \ 

''CR, 

I  J 

AV  _  Is  changed  until 
no 

60 

- 

(v».  ♦  A 

Vno)  ^  toltronct 

Figures  2  and  3  show  the  velocity  versus  axial  chord  distance  plots  obtained  from  the  former 
computer  program  and  the  present  method  for  Identioal  Input  gas  stream  conditions. 


The  computer  program  presently  in  use  at  AFAPL  calculates  convection  heat-transfer 
ooeftlolents  and  adiabatic  wall  temperatures  along  the  blade  profile.  Statlr  temperature  and 
pressure  are  calculated  outside  the  boundary  layer  from  channel  flow  theory*  using  the 
following  equations: 


"s 


Cp  *  /i.  *  Pj^»  and  p  •  characteristic  of  the  fluid*  are  entered  as  tables  for  the 
appropriate  range  of  temperatures  and  pressures. 
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Velocity 
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Hw  calculation  of  convection  heat^tranafar  coefflotonta  and  adiabatic  wall  temperatures 
la  then  carried  out  on  the  following  assumpdona: 


1.  On  the  suction  side*  die  boundary  layer  la  laminar  fr<Mn  the  stagnation  point  to 
the  point  of  minimum  pressure.  From  this  latter  point  on*  the  flow  Is  turbulent. 

2.  On  the  pressure  side*  the  bounckry  layer  Is  turbulent  brom  the  neighborhood  of 
the  stagnation  point  to  the  trailing  edge. 

The  adiabatic  waU  temperatures  are  obtained  from  the  following  equations: 

(lominar)  ^  *"^0  "  V 

(turbulent)  *  Tg  ♦  (T^  "  Ts  * 

On  the  leading  edge  and  In  the  laminar  part  of  the  boundary  layer*  Squire's  method 
(Reference  3)  for  the  calculation  of  heat  transfer  on  a  csdlnder  and  on  a  flat  plate  Is  used. 

In  the  turbulent  part  of  the  boundary  layer*  an  approximation  of  the  von  Karman  formula  for 
heat  transfer  In  turbulent  flows  Is  used. 

«.  •  i  '1 

assuming  the  local  coefHclent  of  skin  brlctlon  for  a  flat  plate  to  be 


where 


(10) 

(II) 


The  fact  that  the  Reynolds  numbers  anticipated  In  the  testing  of  the  blades  will  be  low* 

»  10^*  and  that*  for  the  pressure  profiles  obtained  (Figure  4)*  Is  small  In  an 

extended  neighborhood  of  the  minimum  pressure  point  malces  questionable  the  assumption 
that  the  latter  coincides  with  the  transition  point.  Also*  Equation  11  suggests  tlmt  dM  Is 
rather  Insensitive  to  fluctuations  In  U*  the  latter  being  overslmdowed  by  those  In  x.  A 
comparison  of  Figures  4  through  6  shows  that  the  drastic  changes  In  the  velocity  profiles* 
particularly  on  the  pressure  side*  are  not  reflected  in  d(i)  and  that*  for  Instance*  the 
momentum  thickness  of  the  turbulent  boundary  layer  is  approximately  the  same  on  both 
(nrossure  and  suction  surfaces  beyond  a  surface  distance  of  1.3  Inch*  corresponding 
approximately  to  the  transition  point  on  the  suction  side.  Obviously*  this  behavior  of  6  will 
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be  reflected  in  N^,  since 


^  I 

N,(  vorict  with  9^,  ond ,  therefore  ,  h  (x) 


In  an  attempt  to  evaluate  the  discrepancies  that  might  result  from  the  preceding 
observations!  we  have  adopted  the  theories  of  H.  Schlichting  and  E.  Truokenbrodt (References  4 
and  Si  respectively)  relative  to  the  determination  of  the  transition  point  and  to  the  computation 
of  a  two-dimensional  turbulent  boundary  layer  momentum  thlokness»  to  the  specific 
requirements  of  AFAPL  turbine  blade  hot  testii^  investigation  and  incorporated  them  in  a 
computer  program  described  in  this  report.  These  theories  will  be  described  briefly. 

In  the  new  program,  the  laminar  boundary  layer  is  calculated  using  Scpiire's  method 
which  uses  a  simple  quadrature  to  obtain  6  ,  where 


/u» 

II* 


The  displacement  thickness  S*  and  boundary  layer  thickness  8  are  obtained  from  the 
assumptions 


■f-  *  •  'i’  *  tV 


The  shape  factor  of  the  boundary  layer  velocity  profiles  is  deflned  as 

X  .  -SiL 

W  dx 

The  pressure  decreases  for  X  >0  and  increases  for  X<  0  ,  hence  the  point  of  minimum 
pressure  occurs  for  X  ^  0 

Original  form! 


4*.  -i-_Lx  i  4  =(— -  —  X  - 

5  10  120  '  8  \  315  945 


9072 


A  family  of  neutral  stability  curves  deflned  by  the  shape  factor  X  and  representing 

the  variation  of  8  *  with  the  Rejmolds  number  u8*  has  been  obtained  by  H.  Schlichting 

V 

and  A.  Ulrich  (Reference  4).  The  point  on  these  curves  at  which  the  Reynolds  number 
has  its  smallest  value  is  defined  as  the  limit  of  stability  for  the  laminar  flow  of  interest. 
This  Rg  *  is  called  the  critical  Reynolds  number, 
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It  followii  that  can  be  plotted  aa  a  ftinction  of  X  ,  hence,  of  the  surface  distance 

X.  The  intercept  of  this  curve  with  that  representing  the  variation  of  the  local  Reynolds 
number  with  a  will  3rield  the  separation  point. 

As  may  be  seen  from  Figures  7  and  8  for  small  x's»  R  ^  *  is  small  and  X  is  large.  As  x 
becomes  larger,  this  relationship  is  inverted  and  consequently  Rg^^^ decreases  while  R3  * 
increases. 

The  tormented  pixrfile  of  the  variation  of  R»*  with  x  is  to  be  attributed  to  the  dependence 
dU  ® 

of  X  undergoes  abrupt  changes  along  the  blade  surtece. 

The  calculation  of  6  from  E.  Truckenbrodt's  method  (Reference  5)  is  based  on  the 
energy  Integral  equation 


_!_  A.  '•'»  X '0' 


where  the  energy  thickness  8**  is  a  measure  of  the  kinetic  energy  loss  resulting  from 


friction  and  is  defined  as 


U*  8***  U  (y)  [u*  (I)  -  U*  (y)]dy 


1.12  X  i0~* 

»)* 

rei»raseuts  a  good  approximation  of  the  friction  work  perlormed  in  the  boundary  layer  by  the 
shearing  stresses  r  (H.  Schlichtlng). 

Integration  of  Equation  13  was  performed  by  E.  Truckenbrodt  (Reference  5)  and  yields 


«  C.  ♦  A 


4»sf  ^ 
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which  may  be  rewritten  uainf  the  turbulent  coefficient  of  skin  friction  for  a  flat  plate  as 


represents  the  laminar  portion  of  the  boundary  layer. 

It  may  be  seen  from  Equation  14  that  B  is  highly  sensitive  to  the  variations  in  the 
velocity  U.  Figures  4  through  6  afford  a  comparison  of  the  variations  of  6  and  u  with 
respect  to  x  for  the  suction  and  pressure  sides. 

A  very  good  value  of  the  local  coefRcient  of  skin  friction  could  have  been  obtained  throu^ 
the  relationship 


but  unfortunately  E.  Truckenbrodt's  method  for  the  variatioo  of  the  shape  factor  H  in  the 

range  of  transition  could  not  be  used  because  in  this  range  was  too  small.  Insteadi 

*|c^'  was  obtained  from 

und  this  value  introduced  into  the  von  Karman  formula 

which  takes  into  account  the  variation  of  the  Prandtl  number. 


16 


AFAPL-TR-68-143 


Hie  convection  heat-tranefer  coefficient  for  the  laminar  part  of  the  boundary  layer  is 
calculated  from  S<iuire'8  method*  using  the  universal  function 

obtained  from  intergratlon  of  the  energy  equation  with  <  1.  Squire's  paper  (Reference  3) 
is  referred  to  for  the  details  of  the  procedure.  Tte  new  computer  program*  instead  of 
using  a  table  of  H  versus  A  *  computes  H  ( A)  at  every  point  where  h  is  calculated. 
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SECTION  m 

RESULTS  AND  CONCLUSIONS 


On  the  suction  sldOi  the  transition  point  as  determined  In  the  new  larogram  Is  shifted  by 
0. 37  Inch  downstream.  At  this  point  h^ie  35%  higher  than  the  h  obtained  from  the  former 
program  at  the  same  location  (Figure  9).  On  the  pressure  slde»  Figure  10  shows  that  the  new 
convection  heat-transfer  coefficient  assumes  lower  values  Immediately  downstream  of  the 
transition  point  with  a  maximum  change  of  40%  with  respect  to  the  former  h  at  the  same 
location  and  Increases  at  a  faster  rate  to  assume  a  value  at  the  trailing  edge  representing 
an  Increase  of  36%  over  the  former  h. 

The  average  value  of  h  over  the  blade  remains  essentially  unchanged  when  Truckenbrodt's 
method  Is  used,  so  tnat  the  total  amount  of  coolant  required  for  the  blade  will  not  be  changed 
significantly. 

The  large  discrepancies  In  local  heat-transfer  coefficients  would  necessitate  a 
redistribution  of  the  coolant  flow  formerly  determined  from  the  Allison  program. 

Confirmation  of  the  validity  of  the  Improvements  suggested  in  this  report  will  necessitate 
experimental  data  which  will  become  available  when  the  AFAPL  in-house  testing  program  is 
Initiated. 
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APPENDIX 

COMPUTER  PROGRAM 
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PROGRAM  FUNCTIONS 

LUMIT 

This  is  the  main  program.  It  Is  divided  into  two  basic  parts:  laminar  boundary  layer  and 
turbulent  boundary  layer.  Initially  the  program  reads  in  the  input  data,  calculates  the  initial 
conditions  of  the  laminar  boundary  layer,  iterates  throu^  the  laminar  calculations  until 
transition  Is  reached,  then  proceeds  into  the  turbulent  section  and  Iterates  until  the  end  of  the 
blade  is  reached.  All  integrations  are  performed  in  the  program  and  are  basically  trapezoidal. 

VLOCT  (alternate  entry  CURVE) 

This  routine  reads  In  the  velocity  profile  versus  surface  distance.  It  uses  the  alternate 
entry  to  find  the  velocity  and  the  first  and  second  derivatives  for  a  given  surface  distance. 

It  should  be  noted  that  the  first  and  second  derivatives  are  somewhat  Inaccurate  since  they 
are  based  upon  a  finite  length  rather  than  a  point. 

Uses  subroutine  TLOCK 

DTFRMX 

This  routine  solves  the  universal  function  of  DELTA  for  DELTA  greater  than  1.  It 
iteratively  solves  for  the  unknown  DELTA  given  H  (DELTA). 

Uses  subroutine  AFQUIR. 

See  Chapter  XII  of  Reference  3. 

PRANX 

This  routine  looks  up  in  a  table  the  constant  used  in  the  heat-transfer  equation  near  the 
stagnation  point  as  a  function  of  the  Prandtl  number. 

Uses  subroutine  SRCHX 

PROPEX 

This  routine  calculates  the  viscosity,  density,  Prandtl  munber,  and  specific  heat  of 
air  as  a  function  of  temperature  and  pressure. 

Uses  subroutine  PRCX^OM  and  TLCX)K. 
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TLOOKX 

The  ftinctlon  subprogram  TLOCXC  is  a  general  purpose  routine  to  perform  a  table  look-up 
in  a  two-dimensional  table  (dependent  variable  versus  independent  variable).  It  first  locates 
the  input  independent  variable  in  its  table*  then  takes  the  nearest  'N*  pairs  of  points  and  calls 
subroutine  LAGRNO.  This  program  uses  an  Interpolating  polynominal  of  degree  in 

the  Lagranglan  form*  to  evaluate  the  dependent  variable.  TLOOK  also  has  the  capability  of 
remembering  where  it  found  the  Independent  variable  in  the  table.  Thus*  search  time  is 
saved  when  the  next  time  it  is  called  the  independent  variable  has  changed  only  slightly. 

Use  subroutine  LAORNO. 

LAORNO 

See  description  of  TLOOKX. 

SRCHXX 

This  subprogram  is  a  table  lookup  routine  using  linear  Interpolation. 

PROCOM 

This  routine  calculates  the  thermodynamic  properties  of  air  or  idr-JP4  mixtures.  Given 
temperature  and  ftiel-air  ratio,  it  calculates  speed  of  sound*  ratio  of  specific  heats*  specific 
heat  of  constant  pressure*  gas  constant*  and  nonpressure  biased  entropy  and  enthalpy. 

AFQUm 

This  program  is  a  quadratic  convergence  routine.  It  is  a  routine  having  general 
apidlcation  and  is  used  to  converge  practically  any  function. 

CRITXX 

This  routine  takes  an  input  shape  factor  (ALAM)  and  looks  up  on  a  curve  the  critical 
Reynolds  number. 

Uses  subroutine  SRCHX. 

See  Chapter  XVn  of  Reference  3. 

INPUT  VARIABLES 

NUMB  -  number  of  points  in  velocity  profile  curve 
XS  -  surfoce  distance 

VS  -  surfoce  velocity 

ALENTH  -  characteristic  blade  length 
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TTZERO  -  total  temperature  of  free  stream 
PTZERO  -  total  pressure  of  free  stream 
UCRIT  -  sonic  velocity  of  free  stream 
AF  -  average  specific  heat  ratio 
UIA  -  free  stream  velocity 
DIA  -  leading  edge  diameter 
DX  -  integration  interval 
PRINT  -  print  interval 

OUTPUT  VARIABLES 
Laminar 

X  -  surface  distance 
UX  -  surface  velocity 
T  *  static  temperature  of  surface 
P  -  static  pressure  of  surface 
VISM  -  kinematic  viscosity 
DEL  -  laminar  boundary  layer  thickness 
THETA  -  momentum  loss  thickness 
DELBTR  -  boundary  layer  displacement  thickness 
UNDVFS  -  characteristic  surface  Reynolds  number 
TRMCRT  -  critical  surfMe  Reynolds  number 
ALAM  -  velocity  profile  shape  factor 
H  -  universal  function 
DELT  -  thermal  boundary  layer  thickness 
HL  -  convection  heat-transfer  coefficient 
TND  -  adiabatic  wall  temperature 


Turbulent 


X  -  surface  distance 
UX  -  surface  velocity 
T  -  static  temperature  of  surface 
P  -  static  pressure  of  surface 
VISN  -  kinematic  viscosity 
TAP  -  adiabatic  wall  temperature 
CF2  -  surface  friction  coefficient 
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THETA  -  momentum  lose  thlckneae 
PN  -  Prandtl  number 
RHO  -  denaity 

TUQV  -  Raynolda  number  aaaociated  with  momentum  loaa  thlckneae 
HX  -  turbulent  heat-tranafer  coefficient 

INPUT  FORMAT 

Firat  Garda 

Curve  title  card 

Uaed  to  identify  the  particular  velocity  profile 

Second  Garda 
NUMB 

Number  of  pointa  in  velocity  veraua  aurface  diatance  curve 

Next  NUMB  carda 
XS  VS 

Surface  diatance  and  velocity  carda 

Next  Card 

Caae  title  card 

Uaed  to  identify  the  particular  caae 
Next  Card(a) 

NAMELIST/INPUT/ 

Caae  input;  aee  input  variablea  for  definition. 


(12A6) 


(IS) 


(2F10.0) 


(12A6) 
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lISFTC  LUMIT  M94,XRT 

DIMENSION  ATITLEI12I 

NAMELIST/ INPUT/ALENTHtTT2ERO,PT2EROtUCRIT«AK,UIN,OI At 
AOX.PRINT 
CALL  READ  IN  UMAX  I 
RE ADIS. 5 10 M ATI TLEI I  I. 1-1. 121 
NRITEIAtSllMATlTLEIl  1.1-1.121 

510  F0RMATI12A6I 

511  F0RMATnH1.12A6l 
REA0I5.1NPUTI 
R-OIA/2. 

AINT-0.0 

BINT -0.0 
VISKl-0.0 
URAT 11-0.0 
URATl 2-0.0 
X-0.0 
UX5-0.0 
1-0 

lBlN-1 

ALAM0-T.053 

. . . 

c*****************************  •— — — 

c*************************-***  START  OF  LAMINAR  ••••—•• 
C*******************  —  *******-  .•••—— 

c***-«***. . . 

c 

C  X-'.ERO  CALCULATIONS 
C 

}XB-OX 
*L  AM-7. 05  29 
X«o.O 

CALL  CURVElX.OXBtUX.UPX.UPPXI 

1-1T2ER0 

P-PT2ER0 

CALL  PROPERIT.P.AMU.RHO.PN.CPI 

VISK-AMU/RHO 

VISKO-VISK 

CFL- 1 . 32S/SQRT  lUl MALENTH/ VI SKO I 

OEL-SORTIALAMPVISK/UPXI 

THCTA-0.11T46*OEL 

UELSTk-2.55A*THETA 

CALL  PRANIPN.CONSTII 

ANO-2 . *CONST 1  * SOR T I U1 N*OI A/ V I SK I 

0ELT-2.A01A/AN0 

RAT-OEL/OELT 

H-0.3-0.3*RAT'»0.13333*IRAT*«2.I-0.02142B6*IRAT--A.I 
A«0.005555*IRAT**5.  I 
NL-2.*CP*AMU/IPN*0ELTI*A.62T5 
TAO-T 
HO-H 

OELO-OEL 

C 

C  ENO  OF  X-2ER0  CALCULATl 0««S 
C 

IPS-PRINT/OX'AOX 

UX5S-UX**5. 

NRITEI6.100IX.UX.T 
WRITE  16. 101  IP. VI SK. THETA 
WRITE  16. 1021  DEL. OELSTR.UOOVIS 
WRITE  16. 103 ITRMCRT.ALAM.PN 
WRITEI6.104I0ELTA.0ELT.HL 


HBPK  PME  HANK 
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WRITEI6il0S)T«0 

100  FORNRTIBHO  X>f E 1 3. 7,  SX, 7M  UX>,E13.7t3X«7H  T-«E1S.7) 

101  FORNXTISH  P-i E 1 S.7i SXf 7M  ViSK- . E13.7 t3X •7H  THET«-«E 13 .7) 

102  FORNXTIBH  OE  L- 1  El  3.  7i  3X .  7H0ELSTR- 1  E13.7. 3X«  7HUOOV  IS- lE  1 3  .  • 

103  FORNRTIBH  TRMCRT-f  E 1 3. 7f  3X.  7H  ALRN-.E13.7.3X.7H  PN-iE13.7l 

104  F0RNR7IBH  OELTR-. E13. 7, 3X . 7M  OELT- .  E13.7. 3X.7H  HL-*E13.7I 

103  FORNRTIBH  TA0-fE13.7l 

IP-0 

UX3-UX-P3. 

II  l-l«l 

Rl-I 
X-AI*OX 
IP-IP«1 

14  CALL  CURVEIXfOXfUXfUPX.UPPXI 

IFIUX  .LT.  6.01  CO  TO  11 
UX3-UXPP3. 

T-TT2ER0PI1.-IIAK-1. I/(AK-1. I l•fUX/UCRITI••2.l 

P-PT2ERO*MT/TTiEROI**IAK/IAK-l.ni 

CALL  PROPERIT.P.ANUfRHO.PN.CPI 

VISK-ANU/RHO 

CALL  PRANIPNfCONSTli 

AN0-2.*C0NST1«S0RTIUI N«OIA/VISRI 

0ELT0-2.P0IA/AN0 

lH-1 

0EL02-ALAN0PVISK/UPX 

TERNl-0EL02*UX*«6./l34.*VISKi 

UEINTl-0.0 

0TUM2-0.0 

OELAS-OELO 

0EL2AS-0ELAS*0ELAS 

CO  TO  2 

1  1-141 
Al-1 
X-A1*0X 
1P-1P41 

2  CONTINUE 
UX3S-UX3 
UX$-UX 

CALL  CURVEIX«OXtUX,UPX,UPPXI 
UX3-UX«*3. 

T-TT2ER0*I1.-MAK-1.  I/IAK41.I  l•IUX/UCRITI•P2.  I 
P-PT2EROPIIT/TT2EROIPPIAK/IAX-1.MI 
call  PROPER|T,P,ANUfRHO.PN,CP  I 
VISK-AMU/RHO 
VlSKI-VlSKI«VISK*OX 
TERN1-TERN14MUX3-UX3S  1/2.  ••OX 
C-34.«VI$K/MUX3«UX-UX3S«UXSI/2.I 
0F12A-0EL2AS-C*MUX3«UX3SI/2.  I-OX 
OELA-SORTIOEL2AI 
0EL2E<*TERN1 
0ELE-S0RTI0EL2EI 
OEL-OELE 
OEL2-OEL2E 
THETA-0.11746P0EL 
OELSTR-2.334PTHETA 
UOOVIS-UXPOELSTR/VISK 
ALAN-0EL2*UPX/VISK 
CALL  CRITCLIALAN.TRNCRTI 
IFIUOOVIS  .CE.  TRNCRTI  CO  TO  190 
CO  TO  I30,9«I,IH 
C 

C  USEO  FIRST  TINE  ONLY 
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SO  IH-2 
Xl^X 

CALL  CURVEIXliOXfUXl.UPXltUPMlI 

ZlaPN*l0EL0«UXl«M0t««2./l4.*VISKI 

Z2«Zl«IHO*IUXl'»UXt/2.  t*OX 

ZZ1-0EL0««2.*UX1**6./I34.*VISKI 

ZZZ-ZZl^llUXlAAS.'fUXAAS.t/Z.  ••OX 

M)2-0.1176S«UXM4.«Z2/IPN*ZZ2I 

CALL  OTFMMIH02«OELTA2I 

HG«H02/IOELTA2««2.t 

A0B-HO2/ ( (0. 11765/PNt«UX**4.  I 

CQ0«IH02/I(0.1176S/PN|«UX*«4. II•HG 

TEA«tC«UX«M6«0X 

7EAM0-C/C00 

7EAMB-7EIMO 

TEAMA-B/AOB 

CO  TO  loa 

C 

C  USEO  ALL  BUT  FIBST  TINE 
99  CONTINUE 

TEBNA«TEMA4||UX4UXSI/2.I«0X 
TEANBaTEANB4IIUX««5.4UXS««S.»/2.l*DX 
02HG-  10.  1176S/PN  l•UX••4.•TEItNA/TEANB 
CALL  0TPANI02HG«0ELTA1I 
HC-02HG/I0ELTA1««2. I 

TEMCaTEBNC^I  IUX4UXSI/2.I9I  IHC4HSI/2.  I^OX 
TEANO-TEMB 

02HC1«IO.U76S/PNI«UX«*4.*TEANC/ITERNO*HG» 

CALL  0TPM(02HGl«0ELTA2t 
lOB  H*HG 

0ELTA-0ELTA2 

0ELT-0ELTA90EL 

HL«Z. PCP^ANU/ I PN90ELT 1 *4.62723 
TA0>T4irTZER0>T)«S0RT(PN> 

OXL-OX/ALENTH 

AINTaAINT*!  I  IUX/UINI««S.4|UXS/UIN>«*S.  >/2.  I^OXL 

ClSTR-IICFL/2.  )•SVRTUINTI)••1.2S 

UXS«UX 

HS*H 

IFIIP  .NE.  IPSI  GO  TO  1 
SOO  IP«0 

HRITEIAilOOlX.UX.T 

WRITE  16  1 101 IP. VI SK, THETA 

WRITE  16. IO2IOEL1OELSTR.UOOVIS 

WRITE  I6i103ITRNCRTiALAN. PN 

WRITEI6il04|0ELTA,0ELT.HL 

WRITEI61IOSITAO 

CO  TO  1 

CM***********************************************.*.. 

c*** ••**••••••••• ••••••••••••••  ••••< 

c******************************  ENO  OF  LANINAR  ••••' 

C*********  ••••*•*••••••••••••••  ••••! 

c*************************************************.... 

190  WRITEI6.211IX 

211  F0RNAT(1H2.16HTRANSITI0N  AT  X-,E1S.7| 
WRITEI6.100IX,UX.T 
WRITE  16. lOllP.VISK, THETA 
WRITEI6.102I0EL.0ELSTR.U00VIS 
WRITE  I6.103ITRNCRT.ALAN.PN 
WRITE  16. 10*!0ELTA.0ELT.HL 
WRITE  16. lOSITAO 
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XT>X 

V1SKA>VISK1/XT 

C1STKS>CUTR 

. . . . 

c****************************  ••••••••••• 

c ••••••••••••••••••«•••••••••  START  OF  TUReUL6NT  ••••••••••• 

C* •••••••••••••••••••••••••••  ••••••••••• 

. . . 

IFIIP  .EO.  IPS)  IP>0 
800  ■•■♦I 

AI>1 
UXS«UX 
X>A1  POX 
IP-IPPl 

CALL  CURVEIXiOX,UXfUPX,UPPX ) 

T>TT2ER0*I1.-MAK-1.  )/IAKPl.  ) )  PlUX/UCRITIPPa.) 
P>PTZER0*IIT/TT2ER0)«*(AK/I AK>I.) )l 
CALL  PROPERIT.PtAHU.RHO.PN.CP) 

VISK«ANU/RHO 

TA0>Tp|TTZER0-T)*IPN**I1./3. )) 

BINT-BINTp( ( (UX/UlN)**3.SPIUXS/UlNI**3.SI/2.i*0XL 

CF2>0.0I6/I(UIN*ALENTH/VISK)**0.2SI 

TH0L*IUX/UIN)**I-3. )*(IClSTRSPCF2*BiNTI**0.BI 

THETA>TH0L*ALENTH 

TUOV>THETA*UX/VISK 

TAUOOU>0.012S/(TU0V**0.2S) 

HX«TAUOOU«RHO*U'*UX*|3600./7T8.  )/ 

AII.«5.*S0RT(TAU00U  )*(IPN-1.  )«ALOG(  1.PI5./6.  )*IPN-1.  ) )  i ) 
IFIIP  .EO.  IPSI  CO  TO  900 
823  IFIX  .LT.  XHAXi  CO  TO  803 
nil  STOP 
900  IP>0 

MRITEI6.100fXfUX.T 
WRITE  16, lOllP.VISK. THETA 
HRITE|6,901ICF2>RHO,T«UV 
WRITE I6,902IPN,TA0,HL 


901  F0RMATI8H  CF2«, El  3. 7,5X,  7H  RH0«,E13.7,3X,7H  TU0V*,E13.7I 

902  F0RMATI8H  PN«,E1 3. 7, 3X,  7H  TAO> , E13.7,3X, 7H  HX>,E13.7I 

GO  TO  823 

c*« ••••••«*« •••«•«•••««••«••*«  . . 

c*****************************  ENO  OF  TURBULENT  •••••••••••••••••• 

C«* ••••«••«« •«*•••»•••«••**«••  . . 
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tIBFTC  VLOCT  M94tXR7tOECK 

SUBROUTINE  REAOINIXMAXI 
DIMENSION  XSI200liVS(200l 
DIMENSION  ATITLE(12I 
REAOIS.SIOIUTITLEII  l»I-l«12i 
MR1TEI6,511MATITLE(II«I-I.12i 

510  FORMATI12A6i 

511  FORMATUHl,12A6i 
READIStlOOINUMB 

READIStlOlMXSII  itVSil  l«I-l«NUM8i 

100  FORMATIIBI 

101  FORMATI2F10.0i 
DO  1  I«1,NUNB 
I-l 

1  XSIli«XS(Ii/12. 

XMAX«XS(I  I 
RETURN 

ENTRY  CURVE(X,OX.UEtUEP,UEPPi 

X1«X«0X 

X2-X1«0X 

K«0 

NPT-4 

UE«TLOOK(XtXStVStNUMBtNPTtKtILASTI 

UEl«TLOOK(Xl,XStVStNUMB,NPT.K.UASTI 

UE2«TLOOK(X2tXStVStNUMBfNPTtK,lLASTI 

UEP«((UEl-UEi/(Xl-Xil 

UEPX-  ( (UE2-UE 1  >  /  ( X2-X1 1 1 

UEPP«((UEPX-UEPi/(Xl-Xn 

RETURN 

END 
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tIBFTC  OTFRHX  H9«,XRTtOECK 

SUBROUTINE  OTFRNI ANStOEL I 

DIMENSION  OIBI 

0I2I«0.0 

0I3I-0.0 

AJ-SO. 

TOL  ■0.0001 

OIR-1.01 

TRV-1.5 

1  TRY2«TRY*TRV 

TERM-0. 3«TRY2-0.3«TRY40.1  3)33-10.021  AjaA/TRYJI*' 
A0.00SSSS/ITRY2«TRYI 

CALL  AFOUIRIOiTRYfTERM,ANSfAJ«TOL.OIRtANEH,ICONI 

IFIICON  .EO.  31  CO  TO  S 

IFIICON  .EO.  21  CO  TO  10 

TRYaANEM 

CO  TO  1 

i  NRITEIAilOOIANS 

100  FORMAT  I IHOflOHERROR  IN  OELTA.ANS-.EIS.TI 
10  DEL-TRY 

RETURN 
END 
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tIBFTC  MANX  M94iXKTt06CK 
SUBROUTINE  RRANIR.CI 
OINENSION  RX(9liCX(«l 

DATA  (RX(lltl-lt9l/0.6t0.7.0.Bt0.9*1.0tl>lt7.tl0.tl9./ 
DATA  (CXd  l.l-1.9l/0.4M»0.499»0.921»0.546i0.970t 
A0.992tl*lBil.34.1.94/ 

CALL  SRCMX(R,RX(llt9.0.IL«ATMNI 

C-CX(ILI«ATERN«(CX(IL«1I-CXIILII 

RETURN 

END 
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SIRFTC  PROPEX  M9«tXRTtOECK 

SUBROUTINE  PROPER  (  TtPtAMUtRHOtPNf  CPI 
DIMENSION  HT5(31ltHT6l31l 
C  S  TEMPERATURE  (INOEPENOENT  VARIABLE  TABLE  I 
C  NUMBER  OF  POINTS  31 

DATA  IHTS(KltK-lt31l/«90.t9AO.tA30.«T20.tB10.t900.t990.t  lOeO.t  IITO 

A. , 1260.. 1 3SO. t 1440. t 1930. t 1620.. ITIO. • IBOO. tlVBO. t2160. t 2340. t 2320 

B. .2T00..2Ba0..  3060..  3240..  3420.. 3600.  .3700. .3960.  .4140.. 4320..  4620 

C. / 

C  6  PRANDTLE  NUMBER  VERSUS  TEMPERATURE 
C  IRI 

C  INOEPENOENT  VECTOR  AT  (91 

C  LENGTH  31 

C  NPR  FOR  ORV  AIR 

DATA  (HT6(KI.K-1.31I/.T22..T0B..69T..6B9..683..6B..6B..6B..6B2..6B 
A4.. 686.. 609.. 692.. 696.. 699.. 702.. 706.. 714.. 722.. 726.. 734.. 741.. 749 
8.. 799.. 767.. 703.. 803.. 831.. 863.. 916.. 972/ 

NP>4 

NT-31 

KL-0 

C  T-OECREES  RANKIN 
C  P-PDUNOS/SO.  IN. 

C  ANU-POUNOS6SEC./SO.  FT. 

C  RHO-POUNOS6SEC.  SO. /FT. ••A 
C  CP-FT.SO./(SEC.  SO.  OECREE  RANKINI 
C  TK>BTU/(HR  FT  OECREE  RANKINI 
C  PN-UNITLESS 
C  VIS-POUNOS/HR  FT 
C  T3-0ECREES  KELVIN 
C  RX-BTU/<POUNO  OECREE  RANKINI 
C  CPX-8TU/(P0UN0  OECREE  RANKINI 

CALL  PROCOMI0.O.T.X1.X2.CPX.RX.X3.X4I 
T3-0.9999966T 

VIS-0.00393*T3**1.9/(T3«110.4I 

TK-0.6329*S0RTIT3I*0.00248/I 1.4249.4*10. **(>12. /T3 1 /T3 1 

PN-TL00K(T.HT9.HT6.NT.NP.KL.ILASTI 

CP-CPX*32. 174049*778.26 

ANU-VIS/(  36 00. *32. 1740491 

RH0-(P*144./(RX*778. 26*711/32.174049 

RETURN 

ENO 


34 


AFAPL-TR-ei-143 


SIBFTC  TLOOKX  M94/2.XR7.0ECK 

FUNCTION  TLOOKIF.PT.OT.NTtNP.K,  ILAST) 

C  PT  IS  TABLE  OF  INOEPENOENT  VARIABLES 

C  OT  IS  TABLE  OF  OEPENOENT  VARIABLES 

C  NT  -  SIZE  OF  ABOVE  TABLES 
C  NP  >  NUMBER  OF  POINTS  FOR  INTERPOLATION 
C  K  -  0|  LIMIT  OUTPUT  TO  BOUNDARY  OF  TABLE 

C  K  -  If  EXTRAPOLATE  FOR  VALUES  OUTSIDE  TABLE 

DIMENSION  PTillf  OTIlIt  XTIIOI.  VTIIOI 
IFIILAST  .LE.  01  ILAST  -  I 
I  -  ILAST^l 

IFIPTIll  -  PTINTII  4.15.15 
TABLE  PT  IS  IN  ASCENDING  ORDER 
IFIPTdl  .GE.  P  .AND.  PTIILASTI  .LE.  P)  GO  TO  B 
00  5  I-l.NT 

IFIP  .LE.  PTIIII  GO  TO  B 

5  CONTINUE 
ILAST  ■  NT-1 

6  IFIK  .GE.  II  GO  TO  T 
TLOOK  ■  OTINTI 
RETURN 

T  IL  -  NT-NP^l 

GO  TO  20 

B  IF  1 1  .GT.  II  GO  TO  10 

ILAST  ■  1 

IFIK  .GE.  II  GO  TO  9 
TLOOK  ■  OTIll 
RETURN 

9  IL  -  1 
GO  TO  20 

10  ILAST  -  I-l 
IL  •  I-NP/2 
GO  TO  20 

C  TABLE  PT  IS  IN  DESCENDING  ORDER 

15  IFIPTIll  .LE.  P  .AND.  PTIILASTI  .GE.  PI  GO  TO  B 
00  14  l■l•Nr 

IFIP  .GE.  PTII  II  GO  TO  B 

16  CONTINUE 
GO  TO  6 

20  IFIIL  .LT.  II  IL  •  I 

IFIIL  .GT.  NT-NP«1I  IL  -  NT-NP«1 
00  21  J  ■  l.NP 
L  •  ILPJ-I 
XTIJI  •  PTILI 

21  VTIJI  ■  OTILI 

CALL  LACRNCIP.V.XT.YT.NPI 

TLOOK  •  Y 

RETURN 

END 
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IIBFTC  LACRMG  M94/2,XR7«0ECK 

SUBROUTINE  LAGRNG IX .Vt XTtVT.N) 

THIS  ROUTINE  USES  A  LAGRANGIAN  POLYNOMIAL  BASEO  ON  N  TABULAR 
POINTS  TO  INTERPOLATE  V  AS  A  FUNCTION  OF  X  IN  A  TWO  OINENSIONAL 
TABLE. 

DIMENSION  XTI lit  VTIli 
00  1  I-ltN 

IFIX  .EO.  XTII  It  GO  TO  5 
CONTINUE 
LI  -  1 
L2  -  N-1 

S  -  IXTINI-XTIlll/ABSIXTINt-XTtlll 

6  00  7  I-LltL2 

IFIABSIXTin-XTII^lll  .GT.  .001«ABSIXTM  1 1 1  GO  TO  7 
IFIIX-XTIIII«S  .LE.  0.1  GO  TO  B 
LI  -  I«1 
GO  TO  6 
B  L2  -  I 

GO  TO  9 

7  CONTINUE 
L2  ■  N 

9  V  -  0. 

00  3  I-L1.L2 
Z  -  1. 

00  2  J*LltL2 

IFIJ  .EO.  n  60  TO  2 

Z  ■  Z4IX-XTIJII/IXTIII-XTIJII 

2  CONTINUE 

3  Y  •  Y4ZPYTI I  I 
RETURN 

5  Y  •  YTIII 

RETURN 
ENO 


f 
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tIBFTC  SRCHXX  M94/2,XR7fOECK 

SUBROUTINE  SRCHX I V,VT.Nf KEX, IL tC I 
C 

C  THIS  ROUTINE  LOCATES  V  IN  TABLE  VT 
C  IF  REX  -  Of  LINIT  OUTPUT  TO  TABLE  BOUNDARY 

C  IF  REX  -  1.  EXTRAPOLATE  IF  V  IS  OUTSIDE  TABLE 

OINENSION  VTINI 
IFIVTI2I  .LT.  VTIlll  60  TO  « 

C  TABLE  VT  IS  IN  ASCENDING  ORDER 

00  1  I«ltN 
IFIV-VTini  2i2il 

1  CONTINUE 

IS  IL  -  N-1 

IF  (REX  .EO.  II  CO  TO  3 
C  ■  1. 

RETURN 

2  .L  >  I-l 

IF  1 1  .EO.  11  GO  TO  4 

3  C  -  IV-VTIILII/(VTIIL«1I-VT(ILII 
RETURN 

4  IL  >  1 

IFIREX  .EO.  II  GO  TO  3 
C  >  0. 

RETURN 

C  TABLE  VT  IS  IN  DESCENDING  ORDER 

A  00  7  I-l.N 

IFIV-VTIIII  7,2t2 
7  CONTINUE 

60  TO  IS 
END 
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tIBFTC  PROCOM  M94,XRTtOECK 

SUBROUTINE  PROCOM|FARXfTEXtCSEX,AKEX«CPEX«REXtSEXtHEXI 
IFITEX-300.l2f3f3 

2  NRITEIAtl02l 

102  FORMAT  I lHOf35HPROCOM  INPUT  TEMPERATURE  8EL0N  300.1 
RETURN 

3  IFITEX-4»OO.I9f  9.4 

4  WRITEIA.103I 

103  FORMAT  1 1H0.36HPR0C0M  INPUT  TEMPERATURE  ABOVE  4900.1 
RETURN 

9  IFIFARXI6.7.7 

6  NRITEIA.104I 

104  FORMAT  1 1H0.3BHPROCOM  INPUT  FUEL-AIR  RATIO  BELON  ZERO) 

FARX-0.0 

C  AIR  PATH 

7  CPA  -lllllll.0119940E-29«TEX-1.4924770E-21l*TEX 
1»7.A2197A7E-18I«TEX-1.912B299E-14I*TEX-4.717B376E-12) 

2*TEX«4. 99 194B6E-0BI«TEX-9.1936B7«E-09)*TEX42. 902009  lE-01 

HEA-1 1 1 1 1  M  1. 264442 9E-26«TEX-2.0792922E-22)*TEX 
1«1.2702630E-1BI«TEX-3.029691BE-19)*TEX-1.6794994E-12)*TEX 
242.1B39S26E-OBI*TEX-2.976B440E-09I*TEX42.9020091E-01)*TEX 
3-1.799BBB6E«00 

SEA*42. 902009 1E-01*ALOGI TEX) ( ( 1 1  ( 1.4490767E-26*TEX 

1- 2.42112BBE-22I«TEX«1. 92431 93E-1BI*TEX-3.7B2064BE-19)*TEX 

2- 2.239279OE-12l«TEX«3.2799743E-OBI*TEX-9.1976B79E-09)*TEX 
344.9432300E-02 

1FIFARXI200.200.B 
C  FUEL/AIR  PATH 

B  1FIFARX-.067A23I10.10.9 

9  MR1TEI6.101I 

101  FORMAT  (1H0.A3H1NPUT  FUEL-AIR  RATIO  ABOVE  LIMITS.  IT  HAS  BEEN  RESET 
2T0  0.0676231 
FARX-O. 067623 

10  CPF  ■((  ( (  M7.2678710E-29*TEX-1.333966SE-20)*TEX 
141.0212913E-16I«TEX-4.2091104E-13|6TEX49.96S6793E-10)*TEX 
2-1.3771901E-06I«TEX41.2298630E-03)6TEX47.3B1663BE-02 

HEF«(  ( (HI  (9.0848388E-26*TEX-1.90S0949E-21)*TEX 
141.7021SZ9E-17|6TEX-8.4102208E-14|6TEX42.492169BE-10>*TEX 
2-4.9906332E-07|6TEX46.1293190E-04I*TEX47.3816636E-02> 
3*TEX43.0981930E401 

SEF*47.381663BE-02«AL0G(TEX|4(  I  H  H  1.03B2670E-29*TEX 
1-2.2226118E-21I«TEX42.0429826E-17)«TEX-1.0912776E-13>*TEX 
243.322B928E-10I*TEX-6.BB99909E-07|6TEX41.229B630E-03)*TEX 
346.4B339BE-01 

200  CPEXa|CPA«FARX*CPFI/(1.4FARXI 
HEX«IHEA4FARX6HEFI/(1.«FARXI 
SEX«(SEA4FARX«SEF|/(1.4FARX) 

AMH«2B .97-.9461B64FARX 

REX-1. 9B6379/ANM 

AKEX<PEX/ICPEX-REXI 

CSEX-SeRTIAKEX*REX«TEX*29031.37l 

RETURN 

ENO 
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IIBFTC  AFOUIR  M94,XR7fOECK 

SUBROUTIME  AFQUIR (X,AlNOtOEFENO«ANS.AJ.TOL tOIR tANEM« ICON) 
DIMENSION  X(9I 
C  XI II -NAME  OF  ARRAY  TO  USE 
C  AINO-INOEPENOANT  VARIABLE 
C  DEPEND-  DEPENDANT  VARIABLE 
C  ANS-ANSWER  UPON  MHICH  TO  CONVERGE 
C  AJ-MAX  NUMBER  OF  TRVS 
C  TOL-PERCENT  TOLERANCE  FOR  CONVERGENCE 
C  OIR-OIRECTION  AND  PERCENTAGE  FOR  FIRST  GUESS 
C  ANEW-CALCULATED  VALUE  OF  NEXT  TRV  AT  INOEPENDANT  VARIABLE 
C  ICON-CONTROL  -1  GO  THRU  LOOP  AGAIN 
C  -2  YOU  HAVE  REACHED  THE  ANSWER 

C  -3  COUNTER  HAS  HIT  LIMITS 

C  XI2I-COUNTER  STORAGE 

c  xisi-choose:  method  of  convergence 

C  XI4I-THIRO  DEPEND  VAR 
C  XISI-THIRO  INO  VAR 
C  XI6 I -SECOND  DEPEND  VAR 
C  XI7I-SECONO  INO  VAR 
C  XIBI-FIRST  DEPEND  VAR 
C  XI9I-FIRST  INO  VAR 

C  XI3)  MUST  BE  ZERO  UPON  FIRST  ENTRY  TO  ROUTINE 

Y-O. 

IF(ANS)1.2.1 

1  OEP-OEPENO-ANS 
TOLANS-TOLPANS 
GO  TO  3 

2  OEP-OEPENO 
TOLANS-TOL 

3  IFIABSIOEPI-TOLANSIS.StA 

A  |F(X<2)-'AJI8,B|7 

5  ANEW-AINO 
X(2)-0. 

ICON-2 

RETURN 

6  ANEW-Y 
X(2)-X(2IPI. 

ICON-I 

RETURN 

7  ANEW-Y 
X(2)-0. 

ICON-3 

RETURN 

B  IFIXI3II9,9.12 

C  *•«  FIRST  GUESS  USING  DIR 

9  XI3I-1. 

XIBI-OEP 

XI9I-AIN0 

IFIAINOIlO.ll.lO 

10  Y-OIRPAINO 
GO  TO  6 

11  Y-OIR 
GO  TO  6 

12  IFIXI3l-l.ll3f 13fl6 

C  PPP  LINEAR  GUESS 

13  XI3I-2. 

XI6I-0EP 

XI7I-AIN0 

1FIXIBI-XI6I llAf9f  14 

14  IFIXI9I-XITI|1',9.19 


39 


AFAPL-TR-68-143 


15  A-IXI9I-XI7M/IXISI-XI6II 
V-XI9I-A*XISI 

IFUSSI10.«XI9M-ASSIVII9,9.« 

C  •••  OUAOMATIC  GUESS 

16  XI6I-0EF 
XISI-AINO 

IFIXITI-XISIIlBtlTtlS 
IT  IFIXI6l-XI6inSt9tlS 

IS  IFIXI6l-XI6in9.13.19 

19  IFIXI9l-XI9ll2St20t23 

20  IFIXISI-XI6M21t22«21 

21  XI9I-XI7I 
XIBI-XI6I 
60  TO  13 

22  XI9I-XI7I 
XISI-XI6I 
XI3I-1. 

IFIXI9lll0tlltl0 

23  IFIXISI-XI6ll26,21t24 

26  Fa|XI4l-XI6ll/IXI7l-X(SI I 

A-IXISI>XI6I-F«IXI9I-XISIII/I  IX(9I>XI7I|6|XI9I-XISIM 

a-F-A«IXI9l6XI7ll 

C«XI6|6XI9I«(A«XI7I-FI 

IFIAI262t260«262 

260  IFIBI261,7«261 

261  V—C/B 
GO  TO  37 

262  IFISI267,263t267 

263  IFICI269t266«26S 

266  V-O. 

GO  TO  37 
269  G*-C/A 

IF(6l7t7«266 

266  V-SQRT(GI 
VV«>SOftT(GI 
60  TO  270 

267  lF(CI269t260t269 

26B  V—B/A 

vv«o. 

GO  TO  270 

269  0>6.6A6C/B9«2 
IF(l.-Oll3t29t26 

29  V— B/(2.*A| 

GO  TO  37 

26  E-SOKTIl.-OI 

27  V«(-B/(2.9AII«I1.6EI 
VV«I-B/(2.«AII«(1.-EI 

270  J-6 
0EFMIN-AaS<XI6M 
00  29  I-6faf2 

IF(0EPMIN-ABS<XIIIII29.29,2B 
2B  J-I 

OEFMIN-ABSlXdll 

29  CONTINUE 
K-J61 

IFMXIKI>VI«(X(KI-VVII32f32f30 

30  IFIABSIXIKI-VI-ABS(XIKI-VVII37f37,31 

31  V-VV 

GO  TO  37 

32  IFI^6l33t36,36 

33  JJ-J62 
KKaK62 
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CO  TO  9S 

34  JJmj-2 
KK-K-2 

35  SLOPE-IXIKKI-XIKI 
IFISLOPE*XIJI*IXIKI-VII36.36.37 

36  v»rv 

37  XI9I-XI7I 
X(SI-XI6I 
XI7I«XI9I 
XI6I«XI4I 
CO  TO  6 
ENO 
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SIBFTC  CRITXX  M94(XR7(0ECK 

SUBROUTINE  CRITCLULAM.TERMI 
DIMENSION  ALMTABI13I.TRMTABI13I 

OATAIALMTABII  1 1 1  ■  1 1 13  I /-6.  ,>S.  i -4. ,>3.  t-2.  t>l.  tO., 
Al*t2*t3*t4«t5*f6«/ 

OATAITRMTABII 1 1 1 ■ 1 1 13 1 /0. 1 120. 1 13B. 1 1 75. .250. t3T5. t 
A645 . 1 1 125 . (2000. • 3500. 1  5500  . . 8000. • 10000 ./ 

CALL  SRCHXIALAM.ALMTABI1I.13.0(1L.CI 
TERM-TRMTABI IL l*C«ITRMTABI  1  L^l l-TRMTABI ILI I 
RETURN 
END 
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OAT  A 
FOR 

SAMPLE  PROOLEN 
ONE 

SUCTION  SURFACE 


oso 

0.0 

0.0 

0.0460 

979.01 

0.0920 

1708.05 

0.1380 

1783.37 

0.1B40 

1872.09 

0.2300 

1932.82 

0.2760 

1995.84 

0.3220 

2040.62 

0.3680 

2101.52 

0.4141 

2i56.37 

0.4601 

2244.38 

0.5061 

2311.80 

0.5521 

2366.87 

0.5981 

2422.29 

0.6441 

2472.06 

0.6901 

2533.90 

0.7361 

2611.21 

0.7821 

2653.90 

0.8281 

2708.52 

0.8741 

2789.65 

0.9201 

2841.08 

0.9661 

2875.03 

1.0121 

2  126.40 

1.0581 

2946.46 

1.1041 

2972.77 

1.1501 

2987.53 

1.1962 

2995.23 

1.2422 

2999.11 

1.2882 

1998.92 

1.3342 

2995.76 

1 .3802 

2992.24 

1.4262 

2980.23 

1.4722 

2971.30 

1.5182 

2955.19 

1.5642 

2932.29 

1.6102 

2901.40 

1.6562 

2901.40 

1.7022 

2882.99 

1.7482 

2865.65 

1.7942 

2841.45 

1.8402 

2822.82 

1.8862 

2802.41 

1.9322 

2779.79 

1.9783 

2762.91 

2.0243 

2739.22 

2.0703 

2715.27 

2.1163 

2692.51 

2.1623 

2669.67 

2.2083 

2642.64 

2.2543 

2646.15 

SAMPLE  PROBLEM  NUMBER  ONE 

•INPUT  ALENTH«0.15945,TTZER0«3510.tPTZERO-57.2tUCRIT-2TlA., 
AK-1.26tUlNa900.tOIA«0.013333tPRINT-l.E-2tOX-l.E-S» 
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SUCTION  SUHMCE 
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SAMPLE 

PROBLEM 

NUMBER  ONE 

x« 

D. 

UX* 

0. 

T* 

0.3310000E  DA 

p» 

D.5T2DDDDE  D2 

VISK* 

0.93187TAE-03 

TMETA* 

0.1893AO1E-0A 

OEf 

D.161195AE-D3 

DELSTR* 

0.A839TA6E-0A 

UOOVIS* 

-O.OOOOOOOE-19 

TRMCRT- 

>D.DDDDDDOE-19 

ALAM* 

0.T0S2900E  01 

PN* 

0.77A2500E  DO 

DELTA- 

-D.DDDDDDDE>19 

DELT* 

0.228A993E-03 

HL* 

0.A96819BE  03 

TAD- 

D.351DDDDE  DA 

X- 

D.999DDDDE-D2 

UX* 

0.1T89103E  DA 

T* 

0.333A322E  DA 

R- 

D.AA6123AE  D2 

VISK* 

0.1103231E-02 

TMETA* 

0.31TT282E-0A 

OEL- 

D.27DA991E-D3 

DELSTR* 

0.811AT7BE-0A 

UOOVIS* 

0.1313968E  03 

TRMCRT- 

D.7B59A9SE  D3 

ALAM* 

0.293B866E-00 

PN* 

0.T628I6AE  DO 

DELTA* 

D.1A15952E  D1 

DELT* 

0.3830136E-03 

ML* 

0.2908361E  03 

TAD* 

D.3A87783E  DA 

X* 

D.1998DDDE-D1 

UX* 

0.19A6693E  DA 

T* 

0.330  22ATE  DA 

P* 

D.A255835E  D2 

VISK* 

0.1139101E-02 

TMETA* 

0.5320A30E-0A 

DEL* 

D.A52958AE>D3 

DELSTR* 

0.13988A3E-03 

UOOV 1 S* 

0.2322226E  03 

TRMCRT* 

D.36377A6E  DA 

ALAM* 

0.3068873E  01 

PN* 

0.T61AB53E  DO 

DELTA* 

D.13DA691E  D1 

DELT* 

0.9909T08E-03 

ML* 

C.18T6126E  03 

TAD* 

D.3A83539E  DA 

X* 

D.:997D0DE-D1 

UX* 

0.2090297E  DA 

T* 

0.32T0AT5E  DA 

P* 

D.AD61D36E  02 

VISK* 

0.11T9906E-C2 

TMETA* 

0.6A2T801E-0A 

DEL- 

D.5A72332E>D3 

DELSTR* 

0.16A1660E-03 

UOOVIS* 

0.2918169E  03 

TRHCRT* 

D.S8AAAA6E  DA 

ALAM* 

0.A13T7T8E  01 

PN* 

0.T602210E  DO 

DELTA* 

D.13D533AE  D1 

DELT* 

0.71A3219E-03 

ML* 

0.13AAT96E  03 

TAD* 

D.3A79318E  DA 

X* 

D.3996DDDE-D1 

UX* 

0.22TAa78E  DA 

T* 

0.322629AE  DA 

P* 

D.38D199AE  D2 

VISK* 

0.1 2296 9AE-02 

TMETA* 

0.6B08T16E-0A 

DEL* 

D.5798625E-D3 

DELSTR* 

0.173B9A6E-03 

UOOVIS* 

0.32169T0E  03 

TRMCRT* 

D.7683AS6E  DA 

ALAM* 

0.A873382E  01 

PN* 

0.T38258AE  DO 

DELTA* 

D.13A96A8E  01 

DELT* 

0.7823A0AE-03 

ML* 

0.1A01398E  03 

TAD* 

D.3A733ADE  DA 

X* 

D.A999DDDE>01 

UX* 

0.2A23690E  DA 

T* 

0.31BT962E  DA 

P* 

D.3Sa8D31E  D2 

VISK* 

0.1278979E-02 

TMETA* 

0.T368910E-0A 

DEL* 

D.62735A9E-D3 

DELSTR* 

0.18a2020E-03 

UOOVIS* 

0.3566A65E  03 

TRMCRT* 

D.5AA9121E  DA 

ALAM* 

0.3972960E  01 

PN* 

0.T3613T9E  DO 

DELTA- 

D.1367D97E  D1 

DELT* 

0.85T65A7E-03 

Ml* 

0.12T1T30E  03 

TAD* 

D.3A6799AE  DA 

X* 

D.S99ADDDE-D1 

UX* 

0.25BA312E  DA 

T* 

0.31A386AE  DA 

P* 

D.3353821F  D2 

VISK* 

0.133B901E-02 

TMETA* 

O.T69eA3BE-OA 

DEL* 

D.65A7282E-D3 

DELSTR* 

0.196A138E-03 

UOOVIS* 

0.3T91129E  03 

TRMCRT* 

D.IDDDDDDE  05 

ALAM* 

0.6A89039E  01 

PN* 

0.T33633AE  DO 

DELTA* 

D.139617AE  D1 

DELT* 

0.91A11A5E-03 

ML* 

0.118399TE  03 

TAD* 

D.3A61718E  DA 

X* 

D.6993DDDE-D1 

UX* 

0.2727728E  DA 

T* 

0.3102099E  DA 

P* 

D.31A33S2E  D2 

VISK* 

0.139903 BE- 02 

TMETA* 

0.8105530E-0A 

DEL* 

D.69DD672E-D3 

DELSTR* 

0.2070152E-03 

UOOVIS* 

0.A03613AE  03 

TRMCRT* 

D.IDDDDDDE  DS 

ALAM* 

0.T28169IE  01 

PN* 

0.T3130T0E  DO 

DELTA* 

D.1A12S9DE  D1 

DELT* 

0.97A7819E-03 

ML* 

0.11036T6E  03 

TAD* 

D.3AS9659E  DA 

X* 

0.7992DDDE-DI 

UX* 

0.2869399E  DA 

T* 

0.30385TRE  DA 

P* 

D.2935326E  D2 

VISK* 

0.1A69939E-02 

TMETA* 

0.83128T6E-0A 

DEL* 

D.72A7A68E-D3 

DELSTR* 

0.21TA18BE-03 

UOOV I S* 

0.A23T123E  03 

TRMCRT* 

D.3966627E  DA 

ALAM* 

0.32333IAE  01 

PN* 

0.7AB9302E  DO 
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DELTA- 

TAO- 

0.1A26261E  01 
0.3AAq2A2E  OA 

OELT- 

0.1033678E-02 

HL- 

0.1036072E 

03 

X- 

o.eqqioooE-01 

UX- 

0.2998661E  04 

T- 

0.3030 109 E 

04 

P- 

0.2e05268E  02 

VISK- 

0.19112S3E-02 

THETA- 

0.9204601E-04 

OEL- 

0.78J6371E-03 

OELSTR- 

0.23908SSE-03 

UOOVIS- 

0.4602399E 

03 

TRMCRT- 

0.3313718E  04 

ALAM- 

0.2879812E  01 

PN- 

0.7475999E 

00 

DELTA- 

TAO- 

0.1415q77E  01 
0.344S030E  04 

OELT- 

0.1109612E-02 

HL- 

0.9608484E 

02 

X- 

o.qqqooooE-01 

UX- 

0.299S5S3E  04 

T- 

0.3018067E 

04 

P- 

0.27S16S2E  02 

VISK- 

0.1931149E-02 

THETA- 

0.1015453E-03 

OEL- 

0.864S09SE-03 

OELSTR- 

0. 259346 7E-03 

UDOVIS- 

0.5073881E 

03 

TRMCRT- 

0.9882192E  03 

ALAM- 

.7150400E  00 

PN- 

0.7469943E 

00 

OELTA- 

TAO- 

0.1388053E  01 
0.3443239E  04 

OELT- 

0.1199905E-02 

HL- 

0.8867699E 

02 

TRANSITION  AT  X>  0.10767008^-00 


X- 

0.1076700E-00 

UX- 

0.2998734E  04 

T- 

0.301T021E 

04 

P- 

0.2747037E  02 

VISK- 

0.1532892E-02 

THETA- 

0.1100902E-03 

OEL- 

0.9372568E-03 

OELSTR- 

0.2811703E-03 

UOOVIS- 

0.5500421E 

03 

TRMCRT- 

0.9496133E  03 

ALAM- 

-0.3532841E-00 

PN- 

0.T469456E 

00 

OELTA- 

0.1361217E  01 

OELT- 

0.1275678E-02 

HL- 

0.8340075E 

02 

TAO- 

0.3443083E  04 

X- 

0.1098900E-00 

UX- 

0.2996996E  04 

T- 

0.301T592E 

04 

P- 

0.274999TE  02 

VISK- 

0.1531939E-02 

THETA- 

0.9122335E-04 

CF2- 

0.9145630E-03 

RHO- 

0.7642255E-03 

TUOV- 

0.1784640E 

03 

PN- 

0.7469T22E  00 

TAO- 

0.3464371E  04 

HX- 

0.3184159E 

03 

X- 

0.1198800E-00 

UX- 

0.2977954E  04 

T- 

0.3023830E 

04 

P- 

0.2777209E  02 

VISK- 

0.1521585E-02 

THETA- 

0.1268618E-03 

CF2- 

0.9130137E-03 

RHO- 

0.7703190E-03 

TUOV- 

0.248  2862E 

03 

PN- 

0.7472640E  00 

TAO- 

0.3465007E  04 

HX- 

0.2916595E 

03 

X- 

0.1298700E-00 

UX- 

0.293 55 70E  04 

T- 

0.3037570E 

04 

P- 

0.2838903E  02 

VISK- 

0.1499094E-02 

THETA- 

0.1642179E-03 

CF2- 

0.9096209E-03 

RHO- 

0.7838692E-03 

TUOV- 

0.32  1  5763E 

03 

PN- 

0.74 791 44E  00 

TAO- 

0.3466403E  04 

HX- 

0.2728721E 

03 

X- 

0.1 3986 OOE -00 

UX- 

0.2893645E  04 

T- 

0.3050968E 

04 

P- 

0.2900103E  02 

VISK- 

0.1477577E-02 

THETA- 

0.200 96 74E-03 

CF2- 

0.9063392E-03 

RHO- 

0.7972509E-03 

TUOV- 

0.3935687E 

03 

PN- 

0.7485592E  00 

TAO- 

0.3467759E  04 

HX- 

0.2591173E 

03 

X- 

0.14985 OOE -00 

UX- 

0.2839714E  04 

T- 

0.306T919E 

04 

P- 

0.2979026E  02 

VISK- 

0.1450926E-02 

THETA- 

0.2410173E-03 

AFAPL-TR-68-14S 


CF2- 

0.9022244E-03 

RHO- 

0.01AA22AE-03 

TUOV* 

0.4717127E 

03 

PN« 

0.7494271E  00 

TAO* 

0.3A69A7AE  04 

HX> 

0.2474S87E 

03 

X« 

0.1S98400E-00 

ux> 

0.2786490E  04 

T« 

0.3084336E 

04 

P« 

0.3057077E  02 

V1SK« 

0.142S708E-02 

THETA- 

0. 28216 32E 

-03 

CF2- 

0.89S27e2E-03 

RHO* 

0.8313121E-03 

TUOV- 

0. 551476  7E 

03 

PN« 

0.7S0323SE  00 

TAO- 

0.3471133E  04 

HX- 

0.2377090E 

03 

X« 

0.1698300E-00 

UX« 

0.2732069E  O’ 

T- 

0.3100800E 

04 

P* 

0.3136976E  02 

VISK- 

0.1400985E-02 

THETA- 

0.32S33S1E 

-03 

CF2« 

0.894>3S8AE-fl3 

RHO- 

0.848S095E-03 

TUOV- 

0.6344379E 

03 

PN« 

0.7S12347E  00 

TAD- 

0.3472707E  04 

HX- 

0.2291490E 

03 

X« 

0.1798200E-00 

UX« 

0.2672002E  04 

T- 

0.31 18S9SE 

04 

P« 

0.322S190E  02 

VISK- 

0.1374806E-02 

THETA- 

0.3726592E-03 

CF2- 

0.8901637E-03 

RHO- 

0.8673923E-03 

TUOV- 

0.7242392E 

03 

PN« 

0.7S22297E  00 

TAO- 

0.3474562E  04 

HX- 

0.2211537E 

03 

AFAPL-TR-e8-143 


DATA 

FUR 

SAHPLE  PAOBLEN 
TWO 

PRESSUAE  SURFACE 


052 

0.0 

0.0 

0.020 

650. 

0.0398 

858.61 

0.060 

1220. 

0.0795 

1509.01 

0.1193 

619.57 

0.1591 

597.11 

0.1988 

597.11 

0.2386 

596.67 

0.2786 

597.11 

0.3182 

597.11 

0.3579 

597.11 

0.3977 

602.20 

0.6675 

608.5  3 

0.6772 

626.25 

0.5170 

627.86 

0.5568 

638.86 

0.5965 

657.37 

0.6363 

680.57 

0.6761 

703.61 

0.7158 

722.00 

0.7556 

758.18 

0.7956 

776.07 

0.8352 

818.28 

0.8769 

836.17 

0.9167 

876.76 

0.9565 

896.55 

0.9962 

926.76 

1  .0360 

966. R6 

1.0738 

972.56 

1.1135 

997.11 

1.1533 

1027.38 

1.1931 

1050.58 

1.2328 

1096.16 

1.2726 

1166.06 

1.3126 

1269.38 

1.3522 

1319.13 

1.3919 

1368.27 

1  .6317 

1600.73 

1.6715 

1662.27 

1.5112 

1692.66 

1.5510 

1539.72 

1.5908 

1577.51 

1.6305 

1665.02 

1.6703 

1731.88 

1.7101 

1809.32 

1 . 76  98 

1911.27 

1.7896 

1981.68 

1 .82  96 

2089.07 

1.8692 

2189.85 

1 .9089 

2280.51 

1.9687 

2382.70 

SAMPLE  PROBLEM  NUMBER  TWO 

» INPUT  ALENTH>0.1S9AS,TT2ER0>3S10.tPTZER0*S7.2tUCRIT>271A., 
AK>1.26,UIN>900.,0I A>0.013333t  PR  I MT • 1 .E-2 tOX>l .E-SS 
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PRESSURE  SURFACE 


i' 


l| 


i 
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SAMPLE 


X> 

P> 

OEL« 

tKMCKr. 

OELTA. 

r«D> 


0. 

0.57?OOOOE  0? 
0«lS37866E-03 
-O.OOOOOOOE-19 
-O.OOOOOOOE-19 
0-3510000E  04 


P0O8LEH 


UX« 

VISK- 

OELSTft. 

ALAM« 

CELT. 


»<0H8E«  two 

0. 

0.«3I8774e-03 
0*A613487E-04 
0.7052900E  01 
0>2284S53e-o3 


T« 

tmet*. 

UOOV1S« 

PM« 

Him 


0.3510000E  04 
0«1806377E-04 
'O.OOOOOOOE-19 
0.7742500E  00 
0.A968198E  03 


•  •  (9  I 


*•  O.A940000E-0? 

nti  '  °*^*30857E  02 

THUrar'  °* 2126069e-03 

OeiTA.  0.14971126  01 
0.3494925E  04 


_,;*|*<*  0.1045191 

'ifj?*  0*637805! 

°E<-T.  0.3167091 


X> 

Pm 

CF2. 

PN> 


0.9990000E-< 
0.55580606  ( 
0»8116362E-( 
0*77242816  ( 


X> 

Pm 

CE2> 

PM 


0.1998  0006-01 
0.55684936  02 

0.8113777E_fl3 

0.7725442E  00 


X* 

Pm 

CP2m 

PNm 


0. 29970006 -< 
0.55672036  ( 
0.81140966-0 
0.77252986  0 


X> 

P> 

CE2- 

PN> 


0.39960006-01 
0.5551920E  02 
0-8117885E-03 
0*77235986  00 


X* 

P* 


0*49950006-01 
0*55344366  02 


ux* 

V1SK> 

8H0> 

TAO> 


0*61506996  0 
0*95023866-0 
0*13360086-0 
0.35082906  O' 


ux> 

V1SK> 

8H0> 

TAO> 


0  *  94902886 
0.13379986 
0.35084026 


UX* 

VISK. 

8H0> 

TAO* 


0  *  94917826 
0*13377526 
0.35083886 


UX> 

VISA. 

8N0> 

TAO> 


0.95095246-0’ 
0*13348376-01 
0.35082236  Ox 


UX* 

VISK> 


0.65896146  03 
0.95299236-03 


T« 

TH6TA« 

UbOVISa 

PNa 

Mia 


0.33897446  04 
0*^A972816-04 
0*90379516  02 
0*76536826  00 
0*35435416  03 


Ta 

TH6TAa 
TOOVa 
MX  a 


0*41145646-03 
0*26632736  03 
0*10238146  03 


Ta 

TM6T4a 
TOOVa 
MX  a 


0.34906116  04 
0.48485066-03 
0*30383106  03 
0*95704056  02 


Ta 

TM6TAa 
TOOVa 
MX  a 


0 

0 

0 

0 


*34  90444  6  04 
*50901076-03 
*32029046  03 
*’♦752876  02 


Ta 

TM6T4a 
TOOVa 
MX  a 


0*34884656  04 
0*A6682466-03 
0*30767706  03 
0*10029186  03 


Ta 

TM6T4a 


0*34861956  0  4 
0.42949486-03 


SO 


f 

AFAPL-TR>«8-143 


CFJ-  0.812:236E-03 
PN-  0.77216S9E  00 

*■  O.S994000E-01 
P-  0.5496CiOE  02 
CP2-  0.8i31852E-03 
PN-  0.7717417E  00 

*■  0.6993000E-01 
P-  0.5434495E  02 
CP2«  0.8147424C-03 
PN*  0.7710679E  00 

*■  0.799200CE-01 
P-  0.5377816E  02 
CF2«  0.8161931E-03 
PN-  0.770455 IE  00 

*■  0.8991000E-01 
,^P-  0.5319458E  02 
CE2«  0.817.068E-03 
PN«  0.7698310E  00 

*■  0.9990000E-01 
P-  0.5253218E  02 
CP2»  0.8194486E-03 
PN«  0.7691326E  00 

*■  0.1098900E>00 
$**0*3962E  J2 
CP2.  0.8251241E-03 
PN«  0.76 7005 IE  00 

*■  0.1198800E-00 
P-  0.4911196E  02 
CP2«  0.8288693E-O3 
PN-  0.7659678E  00 

*•  0.1298700E-00 
"•♦^56542E  02 
oi"  0.8333838E-03 
PN-  0.7648301E  00 

*■  0.1398600E-00 

f=^"  02 

CP2«  0.8408131E-03 
PN«  0.7631662E  00 

*■  0.U98500E-00 
rA’  0**»79890E  02 
CP2-  0.8S18649E-03 
PN-  0.7609956E  00 

*■  0«1598400E-00 
°*3^*2«51E  02 
CP2-  0.8671912E-O3 
PN-  0.7578822E  00 


*H0«  0.1331499E-02 
TAO.  0.3508034E  04 

0.7249753E  03 
ifH*  0*’525137E-03 
RH0«  0.1324157E-02 
TAO-  0.3507616E  04 

UX«  0.8203803E  03 
VISK«  0. 9648691 E-03 
PHO-  0.1312371E-02 
tad-  0.3506937E  04 

0*««9e981E  03 
VISK-  0.9717594E-03 
PHD-  0.1301506E-02 
TAO-  0.3506304E  04 

*/!«-■  “•’T56962E  03 

VI SK-  0.9769883E-03 
PHD-  0.i290285E-02 
TAD.  0.3505642E  04 

OX-  0.1055869E  04 
0.9873564E-03 
PHD-  0.1277516E-02 
TAO-  0.3504880E  04 

®*»2«‘>T59E  04 

VI SK-  0.1014995E-02 
PHD-  0.1236960E-02 
TAO-  0.3502390E  04 

„  0*»A08U3E  04 

0*»033550E-02 
PHD-  0.1211049E-02 
TAO-  0.3500756E  04 

?*»*^*309F  04 

®*T05M52E-02 
RHO-  0.1180682E-02 
TAO-  0.3498794E  04 

®**^*T139E  04 

®*1094422E-02 
RMO-  0.1132690E-02 
TAO-  0.34955S3E  04 

OX-  0.2003249E  04 
RhSI  ®*»»53108E-02 
"MO-  0.1065585E-02 
TAO-  0.3490855E  04 

OX-  0.2302305E  04 
®**238359E-02 
RHO-  0. 980301 7E-03 
TAO-  0.3484350E  04 


T- 

THETA- 

TUOV- 

HX- 

T- 

THETA- 

TUOV- 

HX- 

T- 

THETA- 

TUOV- 

HX- 


0.2969809E  03 
0.106 187 IE  03 

0.3481186E  04 
0.3500712E-03 
0.2650542E  03 
0.1197856E  03 

0.3473104E  04 
0.2687294E-03 
0.2284873E  03 
0.13V8232E  03 

0.3465604E  04 
0.2327934E-03 
0.2155784E  03 
0.1545282E  03 

0.3457811E  04 
0.2126662E-03 
0.2119511E  03 
0.1668912E  03 


TmptI"  ®‘’***®82E  04 

®.1979046E-03 
^®®^*  0*2I16372E  03 

MX-  0.1789277E  03 

TmrtI’  ®*®*2®0T4E  04 

0.1360377E-03 
TOOV-  0. 17165741  03 
MX-  0.2224568E  03 

TMerl"  ?*®*®»300E  04 
TMETA-  0.1331502E-03 
TOOV-  0.1814045E  03 
MX*  0.2359600E  03 

tuctI*  ®*3378917E  04 
0*1311976E-03 
TOOV-  0.1920679E  03 
MX-  0.2487954E  03 

THfrl"  ^”*2®5TE  04 
®*“’*315E-03 
TOOV-  0.1909801E  03 
MX-  0.270  1050  E  03 

T-  0.3290000E  04 
THETA-  0.1087852E-03 
TOOV-  0.1889e82E  03 
MX-  0.2921286E  03 

®*3219412E  04 
®*»®»AT62E-03 
rOOV-  0.1886604E  03 
MX-  0.3089611E  03 
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